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Abstract 

We consider the chiral expansion of the octet baryon magnetic moments in heavy baryon 
chiral perturbation theory including all terms which are of order q A . These terms are 
formally of quadratic order in the quark masses. We show that despite the large non- 
analytic quark mass corrections to the Coleman-Glashow relations at order g 3 , including 
all analytic and non-analytic corrections at order q A , which in total are of moderate size, 
allows for a fit to the measured magnetic moments due to the appearance of counter 
terms with free coupling constants of natural size. In this scheme, the AS 0 transition 
moment is predicted to be /x As o = (1-42 ± 0.01) yisr- 
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1 Introduction 


The magnetic moments of the octet baryons have been measured with high precision over 
the last decade. On the theoretical side, SU(3) flavor symmetry was used by Coleman and 
Glashow [IJ to predict seven relations between the eight moments of the p, n, A, XR, X°, S _ , 5° 
and the /iae° transition moment in terms of two parameters. One of these relations is in fact a 
consequence of isospin symmetry alone. In modern language, this was a tree level calculation 
with the lowest order effective chiral meson-baryon Lagrangian of dimension two. Given the 
simplicity of this approach, these relations work remarkably well. Within chiral perturbation 
theory, corrections beyond tree level have been calculated treating the baryons as relativistic 
spin-1/2 fields @@[@]. As it turns out, the large kaon and pion loop contributions leading 
to non-analytic corrections of the type m}J 2 and m q In m q (with m q a generic symbol for 
any of the light quark masses m u ^ s ) tend to worsen the tree level results. With the advent 
of heavy baryon chiral perturbation theory, this problem was reconsidered in ||. In that 
paper, one loop corrections with intermediate octet and decuplct states were considered. It 
turned out, however, that the graphs with intermediate spin-3/2 states did not cancel the 
Goldstone boson loop contributions as it is the case in other observables. It was argued 
that smaller axial couplings are needed to reduce the kaon loops, but no firm conclusion 
was reached. In addition, one relation amongst the magnetic moments free of the axial 
couplings was derived and found to be in good agreement with experiment. Furthermore, 
the baryon magnetic moments were also considered in the large-A'c 1 0 and in a combined 
large-A^c 1 and m s [|7j expansion leading to some interesting relations. In these calculations, 
an extended power counting scheme was used such that graphs with intermediate decuplet 
states start to contribute already at order q 3 . Nevertheless, inspection of fig.2 in [|[] reveals 
that also terms of higher order g 4 (even in the extended counting scheme), which in part 
lead to the non-analytic corrections of the type m q In m q , were included in addition to the 
leading order q 2 and q 3 graphs.^] Recently, a reordering of the chiral expansion based on 
linear flavour symmetry breaking was proposed to overcome the difficulties with the large 
non-analytic corrections^. We will show here that chiral perturbation theory based on 
the conventional counting scheme can indeed describe the magnetic moments with moderate 
corrections beyond next-to-leading order 0(q A ), i.e. they are typically a factor two to three 
smaller than the 0(q 3 ) ones. Our analysis indicates that there is no need to take unphysically 
small axial couplings or include the spin-3/2 states explicitcly. To substantiate this claim, we 
perform a complete calculation at order q A based on the standard power counting, taking into 
account vertices from the dimension two meson-baryon Lagrangian, which are of the type 
(1/mx fixed coefficient), proportional to the scalar second order symmetry breaking and 
also the ones of the form ~ d^cp d„(j), with 0 denoting the Goldstone bosons. The pertinent 
coupling constants related to the latter terms can be estimated by resonance saturation as 
described below. Such contributions have only been considered partly before. If one includes 
the decuplet in the effective field theory and expands the one-loop graphs with intermediate 
spin-3/2 states, one recovers such type of operators (and others which are formally of higher 

#4 As we will show, there are some additional terms at this order related to dimension two insertions with 
fixed coefficients and scalar symmetry breaking. 
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order). We show that there are also important t-channel vector meson contributions which 
previously have not been taken into account. 

The manuscript is organized as follows. The necessary formalism is detailed in section 2, 
with all lengthy formulae relegated to the appendix. We discuss the chiral corrections at 
orders q 2 , q 3 and g 4 , respectively, with particular emphasis on the latter type of contribu¬ 
tions. Section 3 contains the estimation of some low-energy constants. This is based on 
the resonance exchange saturation hypothesis which works well for scattering processes but 
is not yet under control for symmetry breaking terms. Therefore, the low-energy constants 
related to the latter type of operators are left as fit parameters. Our results are presented in 
section 4 and the summary and conclusions are given in section 5. 


2 Formalism 


The starting point of the heavy baryon CHPT ||[|T0[ is an effective Lagrangian formulated 
in terms of the asymptotic fields, here the octet of Goldstone bosons and the ground state 
baryon octet, denoted B. It admits a low energy expansion of the form 


£cff — £ 


M 


r _ r(2) , /*( 4 ) . /-•(I) I /~(2) , r>(3) /'•(h . 

*~-MB M i *■' MB ' MB ' MB ' *~-MB "+ 


( 1 ) 


where the subscript 'M' (' MB') denotes the meson (meson-baryon) sector and the superscript 
'(%)' the chiral dimension. The Goldstone bosons are collected in the familiar SU(3)-valued 
held U = exp {icp/F^}, with the octet decay constant. If not otherwise specified, we set 
F<t> = F k = 1 .21 F n . The baryons are given by a 3 x 3 matrix, transforming under SU(3)/,x 
SU(3 )r as usual matter fields, B —> B' = KBK\ with K(U,L,R) the compensator held 
which is an element of the conserved subgroup SU(3)y. We use the notation of [Til]. The 


(3) 

ellipsis stands for terms not needed here. Note that the terms from C MR are only needed for 
the wave-function renormalization. Beyond leading order, the effective Lagrangian contains 
parameters not hxed by chiral symmetry, the so-called low-energy constants (LECs). In 
principle, these LECs should be pinned down from data or calculated by means of lattice 
gauge theory. We will resort to the resonance saturation method to estimate some of these 
as explained in more detail below. In what follows, we will work to order g 4 in the chiral 
expansion. The pertinent tree and one-loop graphs are shown and catalogued in hg. 1. 
Consequently, the magnetic moment of any baryon takes the form 


.. _ ..(2) (3) (4) 

/i B — /i B + Mb + hs 


( 2 ) 


where the superscript 
contributions. 


(nV denotes the chiral dimension. Let us now dissect the various 


2.1 Magnetic couplings at leading order q 2 

The magnetic couplings of the photon to the baryons start at chiral dimension two and 
the corresponding contribution to the magnetic moments are calculated from the tree graph 
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( 2 ) 

shown in fig. la. These are the terms from C\ IB underlying the Coleman-Glashow analysis 
4“ = ~ 4 SI[F+, B)) - J- bg S-]{F+, B}> (3) 

with S/ the covariant spin-operator, F= —e{v)QF^ v u + uQF^v)) and (.. .) denotes the 
trace in flavor space. Here, Q = diag(2, —1, —1)/3 is the quark charge matrix, u = VU and 
F fll/ the conventional photon field strength tensor. Throughout we work in the isospin limit 
m u = rrid and thus neglect the small A — H° mixing. Notice that we normalize the low-energy 
constants and b® such that we directly get the total, not just the anomalous, magnetic 
moment for the charged particles under consideration. This simplifies the calculation of some 
terms at order g 4 . The Clebsch-Gordan coefficients for the various octet states are collected 
in appendix A, i.e. all octet magnetic moments are given in terms of these two LECs. This 
leads to the Coleman-Glashow relations (we show here the form advocated in [j|]) 

/t-E+ /tp j 2, /tv;— T /in /ip , 

/is- hH” j hH 0 h'n j 2/iyV5H° V^3/in ; (4) 

and the isospin relation 

hs 0 = ^ (^s++ As-) • (5) 


2.2 One—loop contribution at order q 3 

The first loop corrections arise at order q 3 in the chiral counting. They are given entirely 
in terms of the lowest order parameters from Cfj + C^} B . These are F#, the quark masses 
encoded in the meson masses as well as the two axial coupling constants F and D. We use 
here F — 1/2 and D = 3/4. To order g 3 , the one loop calculation leads to corrections of 
the type (we just show the result for the proton since the pertinent coefficients for the other 
baryons are given in appendix A, see also @0), 


4 3) = 


771 

8n 


>9 AL- 2 9 9 -\ / / x 

(D + + -(3 F 2 + D 2 )-^ 


( 6 ) 


One easily convinces oneself that these loop corrections are large, e.g. /A 3) = —2.21,^ 


using the nucleon mass rrijy = 0.94 GeV as the average octet mass. It is important to note 
that for pion (kaon) loops we use the pion (kaon) decay constant. This is legitimate since 
the difference between these two is of order 0{m q ) and thus beyond the accuracy we are 
calculating. As discussed below, the q 3 calculation does not lead to an improvement of 
the simple tree calculation based on eq.([^). Such a behaviour is already known from the 
analogous SU(2) calculation in ref. |TD[ , where e.g. the pion loop contribution ~ M 3 is large 
and cancels 2/3 of the leading term, which is the proton anomalous magnetic moment in 
the chiral limit. It is also obvious from eq.(|6|) that the pion and kaon loops contribute with 

# 5 We give all moments in units of nuclear magnetons (hn) without further specification. 
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similar magnitude. We remind the reader that already Caldi and Pagels [§] derived three 
relations which are valid to this order and independent of the values of the axial coupling 
constants, 

hS+ 2/^E~ j h-H 0 "h P-H - “I - l^n 2/i^ fJ,p , V^P-AE 0 h-E 0 T Fn • (7) 

These are, in fact, in good agreement with the data. The other three relations at this order 
are predictions for the deviations from the Colcman-Glashow relations, the explicit forms 
can e.g. be found in refs. 00. As shown below, these are the ones which pose problems. 
Before one can draw any conclusion about the convergence or failure of chiral perturbation 
theory for the magnetic moments, it is mandatory to get an estimate about the terms beyond 
0(q 3 ). 


2.3 One—loop contribution at order q A 

At order q 4 one has to include one loop graphs with insertions from C^b as we ^ as additional 
tree graphs from jC$ b . Such an analysis has recently been performed for the baryon masses 
and cr-terms ||T0. Note that most of the O(q^) contributions were included in 0, which is 


based on an extended counting scheme. To see this, one has to expand the diagrams with 
intermediate decuplet states of 0 in powers of the decuplet-octet splitting and retain only 
the leading terms. Here, we attempt such a complete calculation for the magnetic moments 
based on the conventional chiral counting with no explicit resonance fields. The fourth order 
contribution to /Xb is given as follows (compare fig. 1) 


a 4 4) = a 4 4>c) + K ( B- d+e+l> + ^ + ^B h+i) + ^b‘ ] 


( 8 ) 


The terms contributing to collect the counter term (tree) contributions with exactly one 
insertion from the dimension four effective Lagrangian. Here, we consider the ones related to 
the explicit breaking of SU(3) due to the large strange quark mass. These have the form|| 


C 


(4) _ 

MB ~ 


^[S",S1 [IF+,B),X + ])- 


ifii 
4 m 


B[S»,S"]B)( X + F+) 


^<B[S^l{{F+,B},X + }> 


( 9 ) 


Here, \ + is th e spurion, x + — clieig(0,0,l), i.e. a factor m s has been pulled out and absorbed 
in the low-energy constants 011 , 2 , 3,4 and /3±. The terms given in eq.(^j) are of chiral dimension 
four since m s = 0(q 2 ) and F pv = 0(q 2 ). The LECs o* and f3 1 will be fitted in what follows 
since the estimation from resonance saturation for such terms is difficult and not yet under 
sufficient control as discussed in detail in ref. [[10. Of course, in general these five terms 
should be written in terms of the full quark matrix, but since m s m,i,7n u: it is legitimate 
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to neglect at this order the pionic contribution. There are two more terms which could 
contribute. These have the form 

d™ = (X+) (x+) (B[S',5"]{F+,B}) , (10) 

where x+ — u^X u ' + U X^ U , X — 2B 0 A4, with M. the diagonal quark mass matrix and 
Bo = mqqm/F^ measures the strength of the spontaneous symmetry violation. We assume 
here B F^. We have already written these two terms in a way which makes obvious that 
the corresponding LECs amount to quark mass renormalizations of b F)F , i.e. 

b%’ F - b°' F + < X+ > b°’ F . ( 11 ) 


Their contribution can therefore be absorbed in the values of the corresponding dimension 
two LECs. We note that the seven terms given in eqs.(p|,pT|) have already been enumerated 
in H (in other linear combinations). 

The next type of graphs are the one-loop diagrams with exactly one insertion from the 
dimension two Lagrangian eq.(||), i.e. they are proportional to the LECs b F,F . For details on 
the calculation of the wave-function renormalization (cf. fig. If), we refer to ref.[|12j]. These 
terms lead to large non-analytic corrections of the form rn q In m q , see also ref. [||]. The various 
contributions for the baryon states are again relegated to appendix A. There are, however, 
three more different types of terms contributing. At this order, we also have to consider 
double-derivative operators at the meson-baryon vertex with the photon hooking on to the 
meson loop, see fig. lg. The corresponding terms of the dimension two Lagrangian can be 
taken from ref. [113 


4L = |S", si i b, (B «„}{«„ B) + ho,11 <B(K, uj, B) ± ) 


( 12 ) 


where (..)-)- means the anti-commutator and commutator, respectively, with the corresponding 
LECs b w and bn, in order. Notice that since u M ~ + 0{4> 2 ), these terms lead to the 

structure d^(f)d v (j) mentioned in the introduction. This contribution to the magnetic moments 
at quadratic order in the quark masses has only been considered partly before. When one 
expands the one loop graphs with intermediate decuplet states, one is lead to operators of 
such type (plus an infinite tower of higher dimension operators). In that case, however, one 
assumes that the corresponding LECs are entirely saturated by these decuplet excitations. 
Our procedure is more general, we will estimate the values of the LECs 69 , 10,11 from resonance 
exchange as detailed in section 3. This includes intermediate baryon excitations as well as 
vector meson exchanges in the t-channel. Furthermore, there are graphs with insertions from 
the part of C^] B which is proportional to 1 /m with fixed coefficients, as depicted in fig. lh+li. 
The corresponding terms of the Lagrangian read 

C-mb = B{B\vD,\vD,B}})-±-{B[D“,\D r ,B}}) 

~ ({ BS „ [£>", {v ■ u, B}]> + IBS „ {v ■ u, p", B]}>) 

~ ({BS, p-, [ii • B]]> + (iS, [« • «, [IT, B]]>) . (13) 
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These terms are most economically constructed by use of the path integral formalism devel¬ 
oped for SU(2) in []TUJ and extended to SU(3) in [fT5| . Of course, they can also be constructed 
with no recourse to the relativistic theory by means of reparametrization invariance |14 . 
These terms have not been calculated before. The explicit expressions of these latter contri¬ 
butions to [ip are also given in appendix A. Finally, there are self-energy graphs (cf. fig. lh) 
with exactly one insertion from the scalar symmetry-breaking terms, 


£ 


( 2 "') 

MB 


= b D (B{ X+ ,B}) + b F (B [ X+ ,B]) 


(14) 


We call the diagrams belonging to this category class (j) graphs. In principle, there is also 
the term ~ bo(BB)(x+)- For our purpose, however, it is sufficient to lump its contribution 
together with the octet mass in the chiral limit. The bo term can only be disentangled if one 
uses additional information like from the pion-nucleon cr-term. The determination of the 
LECs bp, and bp is relegated to the end of the next section. Again, the explicit contributions 
to the fourth order terms in the chiral expansion of the magnetic moments, called /i^ J \ can 
be found in appendix A. We are now in the position to analyse the chiral expansion of the 
magnetic moments. 


3 Estimation of some low-energy constants 


Since we do not take into account explicit resonance degrees of freedom in our calculation, 
these show up indirectly in the values of certain low-energy constants. In particular, it was 
shown in ref. |L5| that in the two-flavor case the LECs of the dimension two Lagrangian can 
be understood in terms of resonance exchange, largely due to the A and t-channel vector 
meson excitations. In the three flavor case, the situation is more difficult since no systematic 
investigation of low-energy constants exists so far. In particular, it is difficult to get a 

m like in eq. (| 


handle on the symmetry breaking terms 


We will therefore not estimate 


these from resonance exchange but rather leave them as free (fit) parameters. On the other 
hand, the LECs 69 , 10 ,n are related to Goldstone boson scattering off baryons. In such a 
situation, resonance exchange can be used easily to estimate the corresponding LECs. We 
follow ref. (T^J and use contributions from the spin-3/2 decuplet, the isovector vector mesons 
and the spin-1/2 Roper octet to fix these. The situation for the LECs b F ’ F is different. As 
shown in eq. (|TT|) , these obtain a piece proportional to the quark masses. Only the quark 
mass independent piece could be gotten from vector meson exchange and one would be left 
with the fourth order pieces as fit parameters. We prefer not to perform this separation and 
work with the total b F ' F as fit parameters. Alternatively, one could constrain the dimension 
two parts from resonance exchange and fit the b e ’ . This is a general phenomenon which 
first appears at order q 4 in the baryon and at order q 6 in the meson sector, respectively. Only 
if one is interested in the behaviour of certain LECs as a function of the quark masses, one 
has to separate the pieces from the various orders. We are not pursuing this issue here. 

Consider now the LECs 69,10,11- We treat the decuplet fields (collectively denoted by A) 
relativistically and only at the last stage let the mass become very large. The pertinent 
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interaction Lagrangian reads 


= f {a - gV)“ 


e^fZJA^e' 


CIJ 


(15) 


with 7 a, 5, c, i, j' SU(3) indices and the coupling constant C can be determined from the decays 
A —> Btc. 0 mi/ (Z) is a Dirac operator, 




(16) 


and Z is an off-shell parameter. With that, the meson-baryon scattering tree graphs with 
an intermediate spin-3/2 decuplet state lead to 


bf = llA, = bn = ^lA, 


with 


r 

/A -T 


m 


+ —^{l + 4Z 2 + 4Z) + 


4Z 


.3(m A — m) ' 3 m\ ' m 4 j 

and ttta = 1-38 GeV the average decuplet mass. In what follows, we use for the coupling C = 
1.2 ... 1.5 |L 6 | and for the off-shell parameter Z = —0.3 ... 0. The larger value corresponds 
to the static A-isobar model and is thus in spirit closest to the study of ref. || where a very 
heavy explicit decuplet contribution was considered. The smaller value stems from the study 
of the P-wave scattering volume CZ 33 |L7| treating the A as a relativistic held (and also from 
pion photoproduction and nucleon polarizabilities). We remind the reader that the pole piece 
due to the decuplet contribution is unambiguous whereas the polynomial pieces depend on 
the off-shell parameters. The latter dependence can always be absorbed in local contact 
terms as it is done here. Similarly, the Roper octet (denoted R) contribution is derived from 
the interaction Lagrangian 


( 17 ) 


(18) 


£rbi, = Ll(/?'."?-(«„. B}) + A(B 7 <* 7 5 [u„, D ) + h.c. 


(19) 


with Dr and Fr axial-vector couplings that reduce in the two-flavor case to Dr+Fr = qa^R 
(with \[r = 0.53 ± 0.04 from the total width of the A*(1440) and ()a = 1-33 from the 
Goldberger-Treiman relation). From the scattering graphs one finds (for more details, see 

0: 

"h Ir , 6 fi = (D\ + Fl) I R , (20) 


4 

3 


b% = --Dll 


R 1 R 1 


Wo — 2 Dr Fl 


with 


Ir = T7T ( -T • (21) 

16 {rriR — m) 

Here, the average Roper octet mass is m R = 1.63 GeV and F R = 0.11, D R = 0.60 as deter¬ 
mined in ref. |12|]. Finally, there is vector meson exchange. We use the vector-meson-pion cou¬ 
pling (V4>4>) based on the antisymmetric tensor held notation of ref. [18|| with \Gy\ = 69MeV 

















from r(p -> 7T7r). The universal vector meson coupling g v is given by gy = g pwn F 2 /My and 
will be used in what follows. For the couplings of the spin-1 fields to the baryons we follow 

ra, 

Cbbv = K d/f ( Bo^ (V pu , B)±) , (22) 

where the two coupling constants Kd/f are constrained by 


K 


D 


Kf = 


9 pNN K-V 


4 y/2 


771 


(23) 


Here, Ky — 6.1 is the tensor-to-vector coupling ratio 


= dfjY,, — d v V, (making use of the equivalence between the vector and 


by reduction to the two flavor case 
of the p and V pu 

tensor field formulations as proven in |TD|). To fix the relative strength of the two tensor 
type couplings, we follow ref.[|2TJ, K D /K F = a M /(l — a M ) with = 3/4. Putting pieces 
together, we arrive at the contribution from the t-channel isovector exchange, 


n V2G v K d Ky l V V2 Gy K f Ky /-. s 

69 -°’ b ^= My = 8m aM ’ = M V ~ 8m 1 ~ QM 


(24) 


We now collect the various contributions and get the following numerical values for these 
LECs for Z = 0 and C = 1.5 (all numbers are in GeV^ 1 ), 


b 9 = + b* + bl = 0.87 - 0.06 + 0.00 = 0.81 , 

b w = &f 0 + bi 0 + b io = 0-44 + 0.02 + 0.50 = 0.95 , 

bn = bf i + 6fx + b\\ = 0.15 + 0.05 + 0.17 = 0.36 , (25) 


which will be used in the next section. We remark that the dimensionless couplings b\ = 2mb l 
are all of order one, i.e. they have ’’natural” size. As expected from the study in ||, the 
decuplet contributes much stronger than the Roper octet. However, whenever nonvanishing, 
the vector meson exchange is also playing an important role. It is important to stress that such 
type of operators were implicitely included in the analysis of [|5[]. This becomes apparent if one 
expands the pertinent Feynman diagrams in powers of A/M w , with A the average decuplet- 
octet mass splitting. This generates (besides others) the contributions b^ lQ , { . What is 
missing in [|] are exactly the vector meson terms shown to be of relevance in Eq.(p5|). We 
remark that this way of treating the LECs induces a spurious dependence on the scale of 
dimensional regularization A. If the LECs were to be fitted from data or resonance saturation 
would be a perfect approximation, all observables would naturally be independent of A. To 
assess this sensitivity, we will later vary A between 0.8 to 1.2 GeV. 

Finally, we have to determine the LECs b F and b F , cf. Eq. (pH|) . Since these appear only 
in the one loop graphs, it is sufficient to perform a fit to the baryon octet masses at order 
q 2 . To be precise, we identify the octet chiral limit mass, including the contribution from 
the term proportional to the LEC bo, with the nucleon mass and fit the two LECs from the 
splittings in the octet. This leads to 

b D = -0.192 GeV" 1 , b F = —0.210 GeV" 1 . (26) 
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Two remarks are in order. First, we note that if one goes to higher orders, the values for 
bo,F are usually different |TT[ The procedure employed here is, however, consistent to 


the order we are working. Second, the identification of the chiral limit mass supplemented 
by the b 0 contribution with the mass of the nucleon is natural since all magnetic moments 
are measured in units of nuclear magnetons, i.e. they refer to 1/rriN- We have also explored 
other options to treat these parameters (the chiral limit mass, bo, bp, bp) but found the one 
described here most satisfactory. 


4 Results 


Consider now the chiral expansion at various orders. A least square tree level fit based on the 
Lagrangian eq.(j3j) leads to b® = 2.39 and 6^ = 0.77. The corresponding magnetic moments 
are given in table 1 in the column labelled 0(q 2 ). The empirical numbers (rounded to two 
digits) are reproduced to less than 20% with the exception of the A and the H’s. This is the 
typical accuracy one expects from such a tree level SU(3) fit. Note, however, that some of 
the Coleman- Glashow relations work much better, for example (/i n + n P + /i£-)//i p = 0 JIJ 
compared to the empirical value of —0.1. 

At the next order, 0(q 3 ), one has large loop corrections (as already mentioned) and thus 
needs to refit the two parameters & 6 ’ . With 

% D = 5.29 , &£ = 2.87 , (27) 


one obtains the values given in the third colum in table 1. Clearly, the description of the 
magnetic moments did not improve, the average deviation from the empirical values is still 
25%. The real problem, however, is that the deviations from the Coleman-Glashow relations 
in most cases go in the wrong direction, i.e. the gap to the experimental values of these ratios 
widens. To be precise, at order q 2 one has the three relations /i p — n%+ = 0, /!=- — = 0 

and /i=o — Hn = 0. At next order, these change to 0.70, 0.83 and 1.67 compared to the 
empirical (mean) values of 0.37, 0.51 and 0.66, respectively. Also, if one calculates the mean 
deviation of the corrected Coleman-Glashow relations to the empirical values, it increases to 
44% at 0(q 3 ) from 22% at leading order q 2 . This is essentially the source of the statement 
that the one-loop chiral perturbation theory calculation fails to describe the baryon magnetic 
moments or that one should decrease the rather large kaon loop contributions as described 
above. 

From the previous considerations it is obvious that one has to go to order q 4 to make 
a definite statement about the convergence and reliability of the chiral expansion of the 
magnetic moments. Some of the non-analytic terms at this order stemming from the insertion 
of the terms in eq.([|) in one loop graphs were considerd in ||. These lead to terms of the 
type m q lnm q , i.e. subleading non-analytic corrections. However, other terms from the 
dimension two Lagrangian can contribute at 0(q 4 ) as it is the case for the baryon masses 
as detailed in [T2|]. These together with the dimension four counter term contributions have 
been calculated here. With seven free parameters, we can of course fit the measured magnetic 
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moments (compare the fourth column in table 1). The LECs take the values 

b% = 3.93 , b% = 3.01 , 

ot\ = —0.45 , «2 = —0.06, a 3 = —0.47, aq = 0.45 , f3\ = —0.74 . (28) 

The values of are changed compared to the 0{q 3 ) fit which indicates that the quark mass 
dependent corrections of order q A are not small. The numerical values of the dimension four 
LECs are of the size expected from naive dimensional analysis since they should be suppressed 
by a factor m s /A x ~ 1/6 compared to the dimension two couplings (with A x ~ 1 GeV the 
scale of chiral symmetry breaking and m s ( A x ) = 175 MeV). Clearly, the predictive power of 
our approach is not very large, only the AE° transition moment is predicted. Our intention 
is, however, more modest. We want to study the convergence of the chiral expansion of 
the magnetic moments. It appears to be under control (we show the results for Z — 0 and 
C = 1.5), 


Ap 

4.31 

(1- 

- 0.51 

± 

0.16) 

Art — 

-2.61 

(1- 

- 0.32 

± 

0.05) 

As+ = 

4.31 

(1- 

-0.64 

± 

0.21) 

As- = 

1 

o 

(1- 

-0.36 

± 

0.05) 

As° = 

1.31 

(1- 

-0.82 

± 

0.32) 

Aa 

-1.31 

(1- 

-0.82 

± 

0.29) 

Ah 0 = 

-2.62 

(1- 

-0.84 

± 

0.32) 

Ah- = 

1 

o 

(1- 

-0.76 

± 

0.15) 

Aae = 

2.28 

(1- 

-0.50 

± 

0.13) 


2.79, 

-1.91, 

2.46, 

-1.16, 

0.65, (29) 

-0.61, 

-1.25, 

-0.65, 

1.42. 


We have set here the scale of dimensional regularization A = 800 MeV. In all cases the 
0(q 4 ) contribution is smaller than the one from 0(q 3 ) by at least a factor of two, in most 
cases even by a factor of three. Like in the case of the baryon masses 0 0 . we find 
sizeable cancellations between the leading and next-to-leading order terms making a precise 
calculation of the 0(q 4 ) terms absolutely necessary. For the E°, the A and the 5° these 
cancellations are so effective that the fourth order contribution appears at first sight fairly 
large. The result for the E° is in agreement with isospin symmetry (as it should) and the 
prediction for the transition moment /tas 0 = 1-42 is within 1.5 standard deviations of the 
recent lattice calculation of [^], /tae° = 1.54 ± 0.09. It is somewhat smaller in magnitude 
than the PDG value, /xaso = ±(1.61 ± 0.08) but within 2.5 standard deviations. To make a 
more realistic prediction for /i s o, one would have to account for effects of the order m u — uid 
and include virtual photons. This goes beyond the scope of this work. 

The various contribution to the fourth order term /i ff are separately given in table 2 (for 
Z = 0 and C = 1.5). Individually, some of these contributions are rather large but in all cases 
sizeable cancellations appear that make the sum relatively small. In particular, we stress that 
the class g diagrams, which were only considered to some extent before, give in some cases 
even the largest contribution. Furthermore, the corrections due to 1/m insertions with fixed 
coefficients are quite sizeable in certain cases like the E + . The terms due to scalar symmetry 
breaking insertions ~ bo,F are in general small, with the exception of the E + and the E _ . 
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This underlines the statement that it is mandatory to take all terms at a given order in the 
chiral expansion. 

It is also instructive to disentangle in the fourth order pieces the various contributions 
from the decuplet, the Roper octet and the vector meson contribution. In all but one cases, 
the decuplet and vector meson contributions are of comparable size and much bigger than 
the one from the Roper octet. Two typical cases are the proton and the neutron (for Z — 0 
and C = 1.5), 


p ■ 

/4 4iA) = 

-0.42 , 

u (4 ’ y) - 

-0.27 , 

i,(4,J?) _ 
r P 

-0.04 

n : 

A) = 

0.23 , 

n (4 ’ y) - 

rn 

0.27 , 

„(4 ,R) - 

r^n 

0.0 


This underlines the previously made statement that the t-channel vector meson excitations 
play an important role in determimning the strength of the insertions ~ &g,10,11 an d that it 
is not sufficient to only take the decuplet to saturate these LECs. 

If one varies the parameters Z and C within the bounds discussed in section 3, the trends 
discussed above are the same. There is some reshuffling between the various fourth order 
contributions, but the sum of them stays approximately the same. Reducing the values for Z 
and C tends to increase the fourth order contribution by some percent. For the proton, the 
relative fourth order contribution varies between 0.16 and 0.20 for the ranges of Z and C given 
above. We note that within these ranges, the prediction for the transition magnetic moment is 
stable within 0.01 p n- For orientation, we give in table 3 the values of the dimension two and 
four LECs for various input values of Z and C. We remark that the relation between certain 
magnetic moments derived in [j5|] does not hold in the complete q 4 calculation presented here. 

Let us now discuss the scale dependence of the results presented so far. As already ex¬ 
plained, our resonance saturation estimate leads to some spurious dependence on A. Typi¬ 
cally, one would vary this scale between the mass of the p and the mass of the A. I11 Fig. 2, 
we show the relative third and fourth order contributions (the second order term is normal¬ 
ized to give one) as a function of A for the proton (which is a representative state). As A 
increases, the fourth order terms become more and more important. However, they stay 
well below the third order contributions. This behaviour is expected since the contribution 
from the chiral logarithms grows as A increases. For example, ln(Af/-<-/().8) = —0.47 where 
as ln(M^/1.2) = —0.88. This increase is only balanced in part by the readjustment of the 
LECs. In particular, the LECs bf and become smaller and thus the leading second order 
contribution decreases which in turn makes the relative third and fourth order contributions 
grow. 

Finally, we remark that there might be potentially large corrections at higher orders. Not 
having performed a complete calculation at 0(q 5 ), the following remarks should be taken 
with caution. We simply note that using the scalar symmetry breaking LECs bn and bp as 
determined from a q 3 |lT| or q 4 pAfl fit leads to very large corrections for the proton from 


the class (j) diagrams which in turn makes the fit values for the a, and (3\ unnaturally large. 
However, there are many other contributions at that order and only a complete calculation 
beyond the one presented here could unravel the true significance of these effects. 
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5 Summary and conclusions 

In this paper, we have considered the chiral expansion of the octet baryon magnetic moments 
at next-to-next-to-leading order 0(q 4 ). Such a study is mandated by the fact that while 
the leading order (g 2 ) relations originally derived by Coleman and Glashow work quite well, 
considering the leading and next-to-leading non-analytic corrections of the type m}J 2 and 
m q In rn q does not improve the description of the magnetic moments and has even some 
times been advocated as a breakdown of baryon chiral perturbation theory for three flavors. 
Our study extends the one of ref. [Q in that we have included all terms at this order .f^] In 
particular, we have shown that the double-derivative meson-baryon couplings generated in 
H by expanding loop graphs with intermediate decuplet states receive important corrections 
from vector meson exchange. In total, we have to deal with seven low-energy constants 
related to symmetry breaking in the presence of the photon held, two LECs related to scalar 
symmetry breaking and three which appear in meson-baryon scattering. The latter can be 
estimated with some confidence from resonance exchange, in our case from decuplet and 
vector meson excitation. This means in particular that we do not include the spin-3/2 fields 
as dynamical degrees of freedom in the effective theory. As stressed in ref.[[L2|, a reliable 
estimation of the LECs related to symmetry breaking is very difficult. We have therefore 
kept the ones related to symmetry breaking in the presence of electromagnetic fields as free 
parameters and are thus able to fit the observed magnetic moments. The two LECs related to 
the scalar sector could be determined from a best fit to the octet mass splittings at order q 2 . 
Our main purpose was to study the convergence of the chiral expansion. We find that with the 
exception of the A, E° and 5°, there appears to be good convergence. In the three exceptional 
cases, the leading and next-to-leading order terms cancel so that the 0(q 4 ) contribution is 
somewhat enlarged. As discussed in detail, there appear large cancellations between the 
various classes of diagrams at order q 4 . In particular, in most cases the contribution from 
the one loop graphs with a double derivative dimension two vertex is substantial and it is 
thus mandatory to include these terms as precisely as possible. The novel ingredient in our 
calculation is the substantial contribution of t-channel vector meson exchange to the LECs 
pertinent to these operators. Only the complete order q 4 calculation performed here can 
be used to explore the convergence of the chiral expansion. The results obtained are stable 
against variations in the input parameters, in particular to the ones related to the decuplet 
sector. We can also predict the AE° transition moment, /iae° = (1-42 ± 0.01)/iat in fair 
agreement with recent lattice gauge theory results. Notice that the theoretical error given 
for Paso refers to the order we have been working and does not account for the uncertainty 
related to higher orders not yet calculated. What still has to be gained is a deeper theoretical 
understanding of the numerical values of the LECs related to symmetry breaking. Work along 
these lines is under way. 


# 6 We stress again that we use the conventional power counting and do not include the decuplet in the 
effective field theory. Our results can be matched onto the ones of || if one expands in powers of (mi — 
win) /M n . How this procedure works in case of the baryon masses is shown in ref. 0- We furthermore 
remark that some 1 /m corrections have not been included in j^] . 
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A Formulae for the magnetic moments 

In this appendix, we collect the explicit formulae for the chiral expansion of the baryon 
magnetic moments to order q 2 , q 3 and q 4 , respectively. The latter are grouped according to 
the various contributions discussed in the main text. 

Order q 2 (compare fig. la) 

/4? = + a F B bl , (A.l) 

with 

Oip = !/3 , of = 1 , «„ = - 2 /3 , «n = 0 , «a = ~1/ 3 , <*a = 0 , 

a° + = 1/3 , a£ + = 1 , o%- = 1/3 , a%- = -1 , ag o = 1/3 , a£ 0 = 0 , (A.2) 

aS- = 1/3 , ag- = —1 , ago = —2/3 , afo = 0 , a^ E o = l/y/3 , a^o = 0 . 


Order q 3 (compare fig. lb) 


with 



M n m K M k m 
8itF 2 + 1 8ttF 2 


(A.3) 


f3; = -(D + F) 2 , (3p = —2/3 (D 2 + 3F 2 ) , /£ = (D + Ff , ft = -(D - F) 2 , 

/5J = 0 , = 2DF ,/3£ + = -2/3 (F 2 + 3F 2 ) , = —(D + F) 2 , 

PI- = 2/3 (F 2 + 3F 2 ) , /f_ = (F — F) 2 , = 0 , /3& = -2FF , (A.4) 

PI- — (D + F) 2 , p%- = 2/3 (F 2 + 3F 2 ) , /3£ 0 = -(F - F) 2 , /3fb = (F + F) 2 , 

/^ s o = -4/V3 DF , = -2/V3DF . 


Order g 4 (compare fig. lc - li) 



mF + Mb' i+e+/> + Mb’’ + Mb'* + ‘ ) 


(A.5) 
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with 


4 4,c) = ai + a 2 + 1/3 a 3 + 1/3 a 4 - 1/3 A , X’ c) = -2/3 a 3 - 2/3 a 4 - 1/3 A , 
/4v’ C) = -8/9a 4 - 1/3 A , /4 4 + } = /4 4 - } = /^S C) = -1/3 A , A® = 0 , (A.6) 

flX = — «2 — 1/3a 3 + 1/3a 4 — 1/3A » = 2/3a 3 — 2/3 a 4 — 1/3 A • 


and 


(4,d+e+f) 

BB 


Ml 


In 




A" 


167t 2 FJ A 


+ &B 


M 2 K 


l6n 2 Fl 


ln 


M; 


K 


X 


+ E ^ 

X=tt, AT, ri 


Ml 

16 n 2 Fl 


In 


M x 

X 


E 




lX=7T, K, T] 


Ml 

8tt 2 FI 



(A.7) 


with 


A 71 " — — h D r) K — —2 h F M — h D 4- h F 

= -2 if , = -6? - = 2 6f , 

= f£ - 6? , «- = 2 6f , & = J? - 6 6 f , 

As 0 = — 2 /a/ 3 6 6 , As 0 = —1/v/l b fi . 


X K — h D _ h F A 77 — n t) K — h D 
°n — °6 °6 > °A — u ) °\ — °6 > 

<£- = E - b% , = 0 , ^ = b° , 

<fo = 6 6 d + b F , (A.8) 


and 


Vp =(D + F) 2 /2 (Jjq - b F ) , = -(F 2 /9 - 2DF + F 2 )b° - (F - F)X , 

r% = -(F - 3F) 2 /18(# + 3& 6 F ) , < = — (F + F)% , 

Vn = (-7/9 D 2 + 2/3FF + F 2 )b F + {D- F)% , r? n = 1/9 (F - 3F) 2 & 6 D , 
hi = -2/3 F 2 &£ , 77* = 1/9 (F 2 + 9F 2 ) b° - 2DFb% , rjl = 2/9 F 2 b° , 
r)l + = 2/9 (F 2 + 6 DF - 6F 2 ) b° - 2F°X , 

E+ = 1/3 ( D 2 + 6FF + F 2 )fe 6 D - (D - F) 2 ^ , r/A = -2/9 F 2 X + 3&£) , 
r/A = 2/9 (F 2 - 6FF - 6F 2 ) b° + 2F% , 

X = 1/3 (F 2 - 6FF + F 2 )^ + (F + F) 2 ^ , rg. = -2/9 D 2 X - O i 
rib = 2/9 (D 2 - 6F 2 ) b» , r/£o = 1/3 (F 2 + F 2 )6^ + 2FF6 6 F , 
nb = -2/9 D 2 b° , ryl- = 1/2 (£> - F) 2 (6 6 D + , 

X- = (-1/9F 2 - 2DF - F 2 ) 6^ + (F> + F) 2 6^ , 

77A = -1/18 (F + 3 F) 2 X - 36f) , 77A — (D — F) 2 b F , 
r/A = (-7/9F 2 - 2/3DF + F 2 )^ - (D + F)X , X = 1/9 (F> + 3F) 2 6^ , 
r/Ao = 2/3/A3 (6FF6g - F 2 6^) , 77^0 = 1/^3 (F 2 - 3F 2 )6f + 2/v^ DFb F , 
r, 7 s „=2/(3V3)D 2 6f. (A.9) 

The wavefunction renormalization coefficients are invariant under SU(2) transformations, 

therefore we give only the values for the different multiplets: 

X = 3/4 (F + F) 2 , X = 5/6F 2 - FF + 3/2F 2 , 0^ = 1/12(71 - 3F) 2 , 
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and 


with 


(pi = 1/3 D 2 + 2 F 2 , <j>* — D 2 + F 2 , p = l/3£> 2 , 

(pl = D 2 , <},f = 1/3 D 2 + 3F 2 , (j)l = 1/3 D 2 , (A.10) 

(PI = 3/4 (D - F) 2 , 0f = 5/6D 2 + DF + 3/2F 2 , (pi = 1/12p + 3F) 2 . 


(A.ll) 


7 p — 2 &io + 2 fen , 7 ^ — &g + 4 6n , 7 ^ — —2 610 — 2 fen , 7 ^ — —2 610 + 2 fen , 

7a = 0 > Ik = —2 610 , 7 s+ = feg + 4 fen , 7 ^+ = 2 bio + 2 fen , 

7s- = — ^9 — 4 fen , 7 *- = 2 & 10 — 2 fen , 7^0 = 0 , 7^0 = 2 b w , (A. 12 ) 

7s- = 2 610 - 2 6n , 72- = -6 9 - 4 6 n , 7 S 0 = -2 &i 0 + 2 6 n , 

7 s° = — 2 6 10 — 2 6n , 7 aso = 4 /\/ 3 feio , 7 as° = 2 /a /3 610 , 


and 


fen = fen + 


8m 


For the graphs with fixed coefficients ~ 1/m, one gets 

3 Ml 


(4 ,/i+i) 
H'B 


= ** In + a* 3M ^ in ^ 

8n 2 F% X +(Jb 8n 2 F 2 X 


K 


Finally, the corrections due to the baryon mass differences read: 

m 


XF = -«b 


2 vr 2 F 2 \ 




nK ™ 

B 2vr 2 Ff 


PP) 


with 


(A.13) 


(A.14) 


(A.15) 


e; = (D + F) 2 (M 2 k b D + (M 2 - Ml) b F ) , 

0f = 1/6 ((3F + DfM 2 + 3 (D - F) 2 M 2 ) 0 D , 

e: = -(D + F)\Ml bn + (M 2 - Ml) b F ) , 0* = {D — F) 2 M 2 b D , 

91 = 0 , 0f = -2 DFMl b D + 1/3 (9F 2 + D 2 ){M 2 K - M 2 ) , 

0£+ = 2/3(3 F 2 M 2 + D 2 M 2 ) b D , 0f + = p + F) 2 {M 2 K b D - (M 2 - M 2 K ) b F ) , 

0£- = -2/3(3 F 2 M 2 + D 2 M 2 ) b D , 0£- = -p - F) 2 {M 2 K b D + (M 2 - M 2 K ) b F ) , 

0-0 = 0, 9 $3 = 2DFMI bn + (A 2 + D 2 ){M 2 K - M 2 ) b F , 

0S- = —P - F) 2 (Ap - {Ml - Ml) b F ) , 

0f- = -1/6 ((3F 2 + D 2 )M* + 3p + FfMl) b D , 

0-o = P - F) 2 (M 2 k b D + {Ml - Ml) b F ) , 0f o = —p — F) 2 M 2 b D , 

0as° = 4/V^3 DF Ml bn , 

0A S o = 2/y/3 DFMl b D + 1 /V3 (3F 2 - D 2 ) {M 2 K - Ml) b F . (A.16) 
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Tables 



0(q 2 ) 

0(q 3 ) 

0(q 4 ) 

Exp. 

P 

2.56 

2.94 

2.79 

2.793 ±0.000 

n 

-1.60 

-2.49 

-1.91 

-1.913 ±0.000 

E+ 

2.56 

2.24 

2.46 

2.458 ±0.010 

E" 

-0.97 

-1.35 

-1.16 

-1.160 ±0.025 

E° 

0.80 

0.44 

0.65 

— 

A 

-0.80 

-0.44 

-0.61 

-0.613 ±0.004 

^0 

-1.60 

-0.82 

-1.25 

-1.250 ±0.014 


-0.97 

-0.52 

-0.65 

-0.651 ±0.003 

AE° 

1.38 

1.65 

1.42 

±1.61 ±0.08 


Table 1: Magnetic moments (in /j,n) calculated to various orders. 


(T) 

/V 

c 

d±e±f 

g 

h±i 

j 

sum 

p 

-0.27 

1.27 

- 0.88 

0.54 

0.01 

0.67 

n 

0.26 

-0.65 

0.47 

- 0.21 

- 0.01 

-0.14 

E+ 

0.25 

1.19 

-0.89 

0.67 

-0.30 

0.91 

E" 

0.25 

-0.56 

0.06 

-0.15 

0.32 

-0.08 

E° 

0.25 

0.31 

-0.42 

0.26 

0.01 

0.42 

A 

-0.15 

-0.39 

0.42 

-0.26 

0.0 

-0.38 


-0.36 

-0.57 

0.72 

-0.53 

-0.08 

-0.84 


0.16 

-0.70 

0.52 

-0.31 

0.08 

-0.25 

AE° 

0.0 

0.50 

-0.44 

0.28 

-0.04 

0.29 


Table 2: Various fourth order contributions to the magnetic moments (in n^). (j) refers to 

diagram (lh) with an insertion ~ &d,f- 


Z 

C 

h° 

h F 

°6 

OL i 

a 2 

OL 3 

(X 4 

Pi 

0 

1.5 

3.93 

3.01 

-0.45 

-0.06 

-0.47 

0.45 

-0.74 

0 

1.2 

3.82 

2.95 

-0.39 

- 0.11 

-0.45 

0.39 

-0.67 

-0.3 

1.5 

3.83 

2.95 

-0.40 

- 0.11 

-0.45 

0.40 

-0.67 

-0.3 

1.2 

3.76 

2.92 

-0.36 

-0.15 

-0.43 

0.36 

-0.62 


Table 3: LECs for various decuplet parameters. 
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Figures 



(a) (b) (c) 

S ' N 



(d) (e) (f) 



(g) (h) (i) 


Figure 1: Contribution to the magnetic moments as explained in the text. The heavy dots 
refer to insertions from c!'m B proportional to the LECs b®, b q, 69 , 10,11 and the crosses refer to 
an insertion with a fixed coefficient. In case of (h), the cross also means an insertion ~ &d,f- 
The box denotes an insertion from C^ B . For the wave function renormalization (f), only one 
representative graph is shown. Also, the crossed partner to (i) is not depicted. 
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Figure 2: Relative contributions to the protons magnetic moment as explained in the text. 
The solid and the dashed-dotted lines refer to the normalized third and fourth order contri¬ 
butions, in order. 


20 





